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Abstract. Since the strange quark has a light mass of order 0{Aqcd), fluctuations of sea ss pairs may 
play a special role in the low-energy dynamics of QCD by inducing significantly different patterns of chiral 
symmetry breaking in the chiral limits TV/ — 2 (rriu = rud ~ 0, rUs physical) and Nf — 3 {rriu = rud ~ 
nis = 0). This effect of vacuum fluctuations of ss pairs is related to the violation of the Zweig rule in the 
scalar sector, described through the two 0{p'^) low-energy constants L4 and of the three-flavour strong 
chiral lagrangian. In the case of signiflcant vacuum fluctuations, three-flavour chiral expansions might 
exhibit a numerical competition between leading- and next-to-leading-order terms according to the chiral 
counting, and chiral extrapolations should be handled with a special care. We investigate the impact of 
the fluctuations of ss pairs on chiral extrapolations in the case of lattice simulations with three dynamical 
flavours in the isospin limit. Information on the size of the vacuum fluctuations can be obtained from the 
dependence of the masses and decay constants of pious and kaons on the light quark masses. Even in the 
case of large fluctuations, corrections due to the finite size of spatial dimensions can be kept under control 
for large enough boxes (L ~ 2.5 fm). 
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In order to achieve a better understanding of nonper- 
turbative features of the strong interaction, it is interest- 
ing to recall the particular mass hierarchy followed by the 
light quarks: 



i^QCD 



H ■ 



(1) 



where Aqcd is the characteristic scale describing the run- 
ning of the QCD effective coupling and Ah ~ 1 GeV the 
mass scale of the bound states not protected by chiral 
symmetry. Therefore, the strange quark may play a spe- 
cial role in the low-energy dynamics of QCD: 

it is light enough to allow for a combined expansion 
of observables in powers of m„, m^, rUs around the Nf = 3 
chiral limit (meaning 3 massless flavours): 



Nf = 3 



mu = rud = iris ^ , 



(2) 



a) it is sufficiently heavy to induce significant changes 
in order parameters from the Nf = 3 chiral limit to the 
Nf = 2 chiral limit (meaning 2 massless flavours): 



Nf = 2 



rriu = md — physical , (3) 



in) it is too light to suppress efficiently loop effects of 
massive ss pairs (contrary to c, 6, t quarks). 

These three arguments suggest that ss sea-pairs may 
play a significant role in chiral dynamics leading to dif- 
ferent behaviours of QCD in Nf = 2 and Nf — 3 chiral 
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limits. Then, chiral order parameters such as the quark 
condensate and the pseudoscalar decay constant: 



SiNf) = - MMuu) , 



F\Nf)^\iTaFl, 

iV f 



(4) 



would have significantly different values in the two chiral 
limits (limjv/ denoting the chiral limit with Nf massless 
fiavours). 

The role of ss-pairs in the structure of QCD vacuum 
is a typical loop effect. Therefore, it should be suppressed 
in the large- iVc limit, and it can be significant only if 
the Zweig rule is badly violated in the vacuum (scalar) 
channel J^'^ = O^^. On general theoretical grounds PP, 
one expects ss sea-pairs to have a paramagnetic effect 
on chiral order parameters. The latter should decrease 
when the strange quark mass is sent to zero : for instance, 
S{2;ms) > S{2;ms — 0), and similarly for F^. This cor- 
responds to 



Z:{2)>E{3), F'^{2)>F'^{3). 



(5) 



However, the size of this paramagnetic suppression is not 
predicted. 

This effect can also be discussed in terms of the Eu- 
clidean QCD Dirac operator, more precisely of its eigen- 
value spectrum (with an appropriate weight over the glu- 
onic configurations) • Chiral order parameters are re- 
lated to the accumulation of the lowest eigenvalues in the 
thermodynamic limit. For instance, the quark condensate 
can be interpreted as the average density of eigenvalues 
around 0. In this language, the Zweig-rule violating ef- 
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feet due to ss pairs eorresponds to multi-point eorrela- 
tions in the density of eigenvalues around [jy. The size 
of the paramagnetic suppression eq. (0) depends on the 
importance of such correlations which can be interpreted 
as fluctuations ^ . 

Thus, it is highly desirable to extract the size of the 
chiral order parameters in Nf = 2 and Nf = 3 limits from 
experiment. Recent data on tttt scattering 0j together with 
older data and numerical solutions of the Roy equations 
allowed us to determine the two-flavour order parameters 
expressed in suitable physical units 6^: 

^(2) = ^^^4^^ =0.81 ±0.07, (6) 

Z(2) = — = 0.89 ± 0.03 . (7) 

A different analysis of the data in ref. with the addi- 
tional input of dispersive estimates for the scalar radius 
of the pion, led to an even larger value of X{2) [7]. In any 
case, X{2) and Z{2) are close to 1, so that corrections 
related to mmrrid (while remains at its physical 
value) have no significant impact on the low-energy be- 
haviour of QCD. In turn, two- flavour Chiral Perturbation 
Theory (^PT) |S] , which consists in an expansion in pow- 
ers of m„ and rrid around the Nf = 2 chiral limit, should 
suffer from no particular problems of convergence. Indeed, 
its two 0{p'^) low-energy constants F^(2) and S{2) are 
dominant in the expansions of the decay constant and 
mass of the pion. 

Unfortunately, two-flavour xPT [5] deals only with dy- 
namical pions in a very limited range of energy. In order to 
include K- and 77-mesons dynamically and extend the en- 
ergy range of interest, one must use three- flavour xPT |S] 
where the expansion in the three light quark masses starts 
around the Nf = 3 vacuum to„ = — rus — 0. From 
the above discussion, large vacuum fluctuations of ss pairs 
should have a dramatic effect on Nf — 3 chiral expansions. 
The leading-order (LO) term, which depends on the 0{p^) 
low-energy constants F'^{3) and would be damped. 

On the other hand, next-to-leading-order (NLO) correc- 
tions could be enhanced, in particular those related to 
Zweig-rule violation in the scalar sector. For instance, the 
Gell-Mann-Oakes-Renner relation would not be saturated 
by its LO term and would receive sizeable numerical con- 
tributions from terms counted as NLO in the chiral count- 
ing. 

We called unstability of the expansion such a numer- 
ical competition between terms of different chiral count- 
ing. A naive argument based on resonance saturation sug- 
gests that higher orders in the chiral expansion should 
be suppressed by powers of (Mt^/Ah)^- However, such an 
argument does not apply to a leading-order contribution 

^ This paramagnetic effect should matter only for observ- 
ables dominated by the infrared end of the Dirac spectrum such 
as the quark condensate and the pseudoscalar decay constant. 
Observables unrelated to chiral symmetry (string tension, vec- 
tor sector) would hardly be affected by this effect and could be 
described accurately through large- A^c techniques. 



proportional to : there is no resonance that could 

saturate the quark condensate. We expect therefore to en- 
counter three- flavour chiral expansions with a good overall 
convergence: 

A = Alo + Anlo + ASA , M < 1 , (8) 

but the numerical balance between the leading order A^o 
and the next-to- leading order Apf^o depends on the im- 
portance of vacuum fluctuations. 

At the level of 0{p'^) Nf — 3 chiral perturbation the- 
ory, the size of the vacuum fluctuations is encoded in the 
low-energy constants (LECs) L4 and Lq whose values re- 
main largely unknown. For a long time, one set them to 
at an arbitrary hadronic scale (typically the 77-mass) as- 
suming that the Zweig rule held in the scalar sector. More 
recent but indirect analyses based on dispersive meth- 
ods [Tnillllll2j suggest values of L4 and Lg which look 
quite modest but are enough to drive the three-flavour or- 
der parameters ^(3) and -^^(3) down to half of their two- 
flavour counterparts -£"(2) and F'^{2). Obvioulsy, these in- 
direct hints of sizeable vacuum fluctuations call for a more 
direct confirmation. 

Unstable Nf ~ 3 chiral expansions {Alo ~ Anlo) 
require a more careful treatment than in two-flavour xPT 
where such unstabilities do not occur. For instance, it 
would be wrong to believe that the chiral expansion oil/ A 
converges nicely ^. This might induce the observed prob- 
lems of convergence in current two-loop computations |13l 
IT^ : the latter treat the fluctuations encoded in L4 and 
Lq as small and are not designed to cope with a large vi- 
olation of the Zweig rule in the scalar sector, leading to 
unstabilities of the chiral series. 

In a previous work we proposed a framework to 
deal with chiral expansions in the case of large fluctua- 
tions, by picking up a subset of observables with (hope- 
fully) good convergence properties and resumming the 
fluctuation terms containing the Zweig-rule violating LECs 
L4 and Lq. This framework includes consistently the al- 
ternatives of large and small vacuum fluctuations. 

Obviously, there is a price to pay for this extension: 
some usual O(p^) relations cannot be exploited anymore, 
because of our ignorance about their convergence. Let us 
comment on a few novelties in our framework: 

— Observables with a good convergence, eq. 0, form 
a linear space, which we identify with connected QCD 
correlators away from kinematic singularities. This choice 
promotes and F'^Mj, {P = tt, K, 77) : LO and NLO 
may compete, but there should be only a tiny contri- 
bution form NNLO and higher. On the contrary, the 
chiral expansion of Mp (ratio of the former quantities) 
may exhibit a bad convergence. 

In principle, two-loop computations could provide a 
check of these assumptions. However, it is difficult to 
exploit available analyses jl3lll4j for two reasons: i) 
in the algebraic expressions, low-energy constants are 
traded for physical quantities assuming that the chiral 

^ This would be equivalent to claim that l/(f + a;) ~ 1 — x 
is a reasonable approximation for x — 0(f). 
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expansions of the latter are dominated numerically by 
the LO contribution, ii) the numerical results rely on 
a specific model of resonance saturation for the 0{p^) 
LECs ^. The dependence on these two assumptions 
should be assessed before any definite conclusion can 
be drawn from two- loop computations ^. Obviously, if 
eq. © were not followed by Fp and FpMp, but rather 
by other combinations of Fp and Mp, some of our 
conclusions might be modified, but this task is beyond 
the scope of the present paper. 

The three-fiavour quark condensate and pseudoscalar 
decay constant expressed in physical units: 



X(3) 



Fue 



Z(3) 



F2 



(9) 



are free parameters. Constraints come from the vac- 
uum stability and the paramagnetic inequalities |SJ: 



< X(3) < X(2) , 



< Z(3) < Z(2) , (10) 



where the values of X{2) and Z{2) have been deter- 
mined from experiment, see eqs. QUI- 
The quark mass ratio r = m,s/m {m = m„ = nrid) 
cannot be fixed from Mj^/M^ since we do not control 
the convergence of its three-fiavour chiral expansion, r 
becomes a free parameter which can vary in the range: 



n = 2 



FkMk 
F^M^ 



p2 i\/r2 

(11) 

The chiral expansions of F^M^ and F'^M]^ lead to 
a correlation between r and X(2) (J. Experimentally, 
the analysis of tttt scattering phase shifts 0| shows that 
r > 14 (95% CL), slightly favouring values between 20 
and 25 [S]. 

The agreement of the pseudoscalar spectrum with the 
Gell-Mann-Okubo formula requires a fine tuning of Lj. 
Let us remark though that this fine tuning exists even 
in the case of a dominant three-flavour quark conden- 
sate and small vacuum fluctuations ,15, . 
One cannot determine LECs or combinations of LECs 
through ratios of observables. For instance, one should 
not use Fk/Ft^ to determine L5 because we do not 
know if the chiral expansion of Fk/Fj^ converges at 
aU. 



^ For instance, in the resonance models used in refs. |13I14) . 
SU{3) breaking in quark masses is taken into account through 
a single constant in the vector sector (fx), and it is neglected 
in the scalar sector {dm = 0). 

* Keeping in mind these issues, it remains an interesting ex- 
ercise to study the convergence of two-loop computations for 
masses and decay constants. In table 2 in ref. |14| . four sets are 
considered which yield rather different results for the overall 
convergence of Fp, Mp, Fp and FpMp {P = tt, K). Set B ex- 
hibits small NNLO terms for Fp, Mp and FpMp, but not for 
In Fit 10, Fp, Fj,, Fj.Ml follow eq. ©, but not M|,. For 
Set A and C, Fp and Mp converge, whereas Fp and FpMp 
suffer from large NNLO corrections. 



There are some prospects of probing experimentally 
the three-flavour sector and in particular the size of vac- 
uum fluctuations through ttK scattering ^21- Unfortu- 
nately, the current data are not precise enough to draw 
any definite conclusion. On the other hand, recent progress 
achieved by lattice unquenched simulations JH] makes 
them an interesting field to investigate the size of vac- 
uum fluctuations. For instance, lattice practitionners can 
vary very easily the value of the quark masses. 

A comment is in order at this point. From large- iVc 
considerations, it is often assumed that quark-loop effects 
are not significant, so that simulations with only two dy- 
namical flavours, or even none (quenched case), could be 
good approximations to real QCD. We are precisely ques- 
tionning this assumption in the case of observables re- 
lated to chiral symmetry breaking. Thus, it is mandatory 
to perform unquenched simulations with three dynamical 
flavours in order to probe strange sea-quark effects on the 
pattern of chiral symmetry breaking. 

One could think of investigating directly the correla- 
tions of Dirac eigenvalues, or the ms-dependence of the 
quark condensate and pseudoscalar decay constant. How- 
ever, this might prove rather challenging tasks since they 
require simulations with three dynamical very light quarks. 
In this paper, we investigate another way of probing the 
size of vacuum fluctuations through the spectrum of the 
theory. Indeed, this less immediate approach is easier to 
follow with current lattice simulations. We believe that 
this exercise might also be useful to illustrate some sub- 
tleties arising in three-flavour chiral extrapolations when 
vacuum fluctuations of ss pairs are not negligible. 

We are not going to address the issue of discretisa- 
tion, which depends on the specific implementation of the 
lattice action. We focus on chiral extrapolations with po- 
tentially large vacuum fluctuations and on the impact of 
finite-volume corrections. For deflniteness, we work in the 
isospin limit to,; — nid — m and we consider a lattice 
simulation with (2-1-1) flavours : two fiavours are set to 
a common mass m, whereas the third one is kept at the 
same mass as the physical strange quark to^. Each quan- 
tity X observed in the physical situation (m, m, mj) has 
a lattice counterpart X for (m, m, mg). 

We introduce the notation: 



m 



m 



X(3) 



2mZ'(3) 

r(3) = 



Z{i) 



FH3) 



F2 ' 
X{3) 2mBo 



Z{3) 



Mi 



(12) 
(13) 
(14) 



where ^(3) and F{3) are (the absolute value of) the quark 
condensate and the pion decay constant in the Nf = 3 
chiral hmit, and Bq = S{3)/F{if. We take the follow- 
ing values for the masses and decay constants: = 92.4 
MeV, Fk/F^ ^ 1.22, = 139.6 MeV, Mk = 493.7 
MeV, A/„ = 547 MeV. 
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1 Vacuum fluctuations at infinite volume 



1.1 Bare expansions of masses and decay constants 



The resulting bare expansions for the decay constants 
(third step of our procedure) can be expressed in the fol- 
lowing way: 



We start by considering the impact of vacuum fluctuations 
on the spectrum in the hmit L — > oo. Since large vacuum 
fluctuations are allowed, the problems highlighted in the 
introduction might arise. Therefore, we have to define the 
appropriate observables to consider and the treatment of 
their chiral expansion. We follow the procedure advocated 
in |15| , reexpressing it in an equivalent way convenient for 
our purposes: 

1. Consider a subset of observables that are assumed to 
have a good overall convergence - we call them "good 
observables" . They must form a linear space, which 
we choose to be that of connected QCD correlators (of 
vector/axial currents and their divergences) as func- 
tions of external momenta, away from any kinematic 
singularities. This rule selects in particular Fp and 
F^APp. 

2. Take each observable and write its NLO chiral expan- 
sion. 

3. In theses formulae, reexpand the physical quantities 
(masses, decay constants. . . ) in powers of quark masses 
wherever the resulting dependence on the latter is poly- 
nomial. Keep the physical quantities only in the nonan- 
alytic terms (unitarity cuts, logarithmic divergences. . . ). 
The result is called "bare expansion" . 

4. In the bare expansions, reexpress 0{p^) and 0{p*) 
LECs in terms of ^(3), Z{3) and r using exact Ward 
identities for pseudoscalar masses and decay constants. 

The first step is straightforward : one has to compute 
the NLO chiral expansion of the masses and decay con- 
stants of Goldstone bosons in an infinite volume with the 
quark masses m„ = rrid = m and rrig. The expressions 
can be obtained directly from ref. [H]. Then, we must re- 
expand the physical quantities occuring in these formulae 
in terms of the quark masses wherever the dependence is 
analytic. In the case of Fp and FpMp, the only issue lies 
in tadpole contributions such as: 



^log^. 



(15) 



To obtain the bare expansions associated with Fp and 
FpMp, we keep the physical masses only in the (nonana- 
lytic) logarithm but we expand the front factor in powers 
of quark masses. Hence, at the chiral order we are working, 
the contribution from the pion tadpole has the following 
bare expansion: 



2mBo . Ml 
log ^ . 



167r2 



(16) 



In the bare expansion, the polynomial structure in the 
quark masses is explictly displayed, so that we can keep 
track of the relative contribution of the LO, NLO, NNLO. . . 
terms. 



F2 

p2 



Z(3) 

F2 327r2 
-Hi 



M2 



qrY{3) 



AP 
41ogH^ 



2 K 



(17) 



F3 



(18) 



TT 

M^arr(3) 



2]l: 



i?2 

V 
p2 



3 

+ 2 

p2 



= zm 

AP 



3 Ml 

2i°s7^ 



1 ) log 



l\Ll 



1 1 log ^ 



(19) 



2 K 



Ml 1 , , /I 



8 /2 
3 [q 



I] LI 



p2 

p2 ' 



where ep are NNLO remainders of 0{m^) (rnq denotes 
either nis or m). We have divided by the physical value of 
Fl in order to deal with dimensionless quantities. 

In a similar way, we obtain the bare expansions of the 
masses: 



FX 
F2MI 



qrl X(3) 



(20) 



M^ 
Ml 



1 



16 



2 a 



F3 327r2 



p2M^ ~ 



qr[Yi3)Y 



AP 



3 log- 

1 Ml, 
log^ - 



1 log ■ 



Fm^ 



FlMl 

Fmi 



qr 
~2 



AP 

-^^-\Ym 




(21) 



Fl 



-8,1 

.9 



1 U£ 
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Ml 



1 



F2 32^2 



qr[Y{i)f 



3 M 2 



+ ■ 



3 (I 



1 log 



Ml 



5 

18 



1 lot 



M2 



ALi = LKfi) — Li{fi) combine the (renormalized and 
quark-mass independent) constants -/j4,5.6,8 f^nd chiral log- 
arithms so that they are independent of the renormalisa- 
tion scale /i: 



327r2L4(M) 



F2M2 



FX 




/ M2 



Ml 



2\Ll 



16 

y 



1 a 



32^2 
1 



qr[Y{3)]' 



4 /2 
9 [q 



1 log- 



i^2 32Tr- 



:qr[Y{3)f 



log^ 

^2 



3 





-I 


(^0 









32 



1 \ M2 

--ij _|:gr[r(3r[3Lr-HL^] 



M2 



where dp are NNLO remainders of 0{rn^). We have di- 
vided by the physical value of F^M^ in order to deal with 
dimensionless quantities. 



1.2 Expression of O(p^) LECs 

We must know perform the fourth step of our procedure. 
As shown in refs. [T1ll()lllllll5j . one can reexpress the 0{p*) 
LECs U, L5, Le, Ls in terms of [r, X{3), Z{3)] (and NNLO 
remainders) using the four chiral expansions of Fp and 
FpMp {P = TT, K) in the physical case. These quanti- 
ties, related to two-point functions of axial/vector currents 
and of their divergences at vanishing momentum transfer, 
are expected to have small NNLO remainders. One ob- 
tains [TS] : 



1 M2. 
-log^ 
8 /i"^ 



(22) 



8(r- l)(r + 2) 
Ml 



(27) 



A/2 

-2 

K 



(4r + l)log^ + (2r+l)log^ 



, Ml M^ 
log^ + 21og^ 



(28) 



r-1) 



3 log 



M2 



32^2Xg(^) = L 
16 



, Ml 2^ Af2 



3 log 



(29) 



16(r + 2)(r-l) 



Ml 
Ml 



M2 
log TTT 



M 



K 



32^2^g(^) = _L 



, Ml 2^ Af2' 



(30) 



1 



16(r- 1) 



01 Ml M2 
3 log — -I- log — f 

Ar2 M2, 



The right hand-side of eqs. involves the r- 

dependcnt functions: 



e{r) 
i]{r) 



To — r 
>_l 

V2-1' 

2 



r2 



FkMk 
F^M^ 



El 
-l\Fl 



1 



1-36, (31) 
(32) 



whereas d, d' and e, e' are combinations of NNLO remain- 
ders associated with the chiral expansions of tt, K masses 
and decay constants respectively. These remainders should 
be small, and we are going to neglect them for all numer- 
ical estimates in the following. 

We stress that eqs. (|23|) - (|26|l are nothing more than 
a convenient reexpression of the Nf = 3 chiral expan- 
sions of F^,Fl,F^MlFlMl. The latter can be easily 
recovered through Hnear combinations of eqs. 
For instance, the chiral expansion of F^ [F^M^] is ob- 
tained when one combines eqs. (|23|I - H24(I [eqs. (|25|l - H26(l ] 
to eliminate ri{r) [e{r)]. 



Y{3)ALi 
Y(3)AL5 
Y\3)ALg 
Y^{3)AL. 



F' 



8 — 



8(r + 2) Ar2 
1 

16(r + 2) Af2 
1 F2 

T6Ay[^« 



[1 - r7(r) - Z(3) - e] , (23) 
e'] , (24) 
■[l~6(r)-X(3)-rf],(25) 
d'] . (26) 



1.3 Numerical results 

To exploit the previous formulae, we have to fix the values 
of r, X{3), Z{3). Since these parameters are only weakly 
constrained by experimental data, we will vary them in the 
ranges r = 20, 30, Z{3) = 0.4, 0.8 and X{3) = 0.2, 0.4, 0.8. 
In each "physical" situation, we study how the masses 
and decay constants of the simulated tt and K vary with 
q, neglecting all NNLO remainders for the moment. 
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A slight computational difficulty arises: in eqs. p7|l - 
H22(l . the logarithmic NLO corrections involve the values 
of the simulated masses Mp. We could solve iteratively the 
system of equations to determine these masses. However, 
two iterations turn out to achieve a sufficient numerical 
accuracy. This corresponds to an easier procedure: i) com- 
pute Mp defined from eqs. (Ejl-ldSll with Mp replaced by 
Mp on their right-hand side, ii) consider eqs. (I17ll - (|22ll 
again, with Mp replaced by Mp on their right-hand side, 
in) the resulting values for the pseudoscalar masses are 
equal to AIfy up to a tiny error. 

In figs. H] and El we have plotted F|Af|, and M|, 
[P = n,K) as functions of q ~ ih/ms in an infinite vol- 
ume, neglecting all NNLO remainders. For each observ- 
able, the first column corresponds to Z{i) = 0.8, the sec- 
ond one to Z(3) = 0.4. On each plot, the curves corre- 
spond to different values of X{2,) [full: 0.8, long-dashed: 
0.4, dashed: 0.2] and of r [thick: r = 30, thin: r = 20]. 

A few comments are in order; 

— Our choice of normalisation imposes that all the curves 
intersect at the physical point m — m where: 



2 Finite-volume effects 
2.1 NLO chiral expansions 

The lattice simulations are performed in a finite spatial 
box, whereas time is sent to infinity to single out the 
ground state. For sufficiently large boxes {L 1/Ah), 
the Goldstone modes remain the only relevant degrees of 
freedom, whose interactions can be described through a 
low-energy effective theory [T^j. In the case of periodic 
boundary conditions, this effective theory is identical to 
xPT, with the same values of the LECs as in an infinite 
volume. In addition, if the box size is large enough com- 
pared to the inverse Compton length of the pion |17II18I 
I19| . the so-called p-expansion is valid and the only differ- 
ence from the infinite- volume case shows up in the propa- 
gators of the Goldstone modes. This affects only the tad- 
pole logarithms in the formulae of the previous section: 



1 
r 



(33) 



If Y(3) > 1, vacuum instability may occur : the pion 
mass becomes negative for "large" masses [q = 0(1)]. 
This situation would occur if vacuum fiuctuations had 
a more significant impact on the decay constant than 
on the condensate : the first would decrease more quickly 
than the latter from the physical case to the three- 
flavour chiral limit. 

Actually, this situation seems unlikely. A nonvanish- 
ing is equivalent to the spontaneous breakdown of 
chiral symmetry. Other chiral order parameters (like 
S) may or may not vanish depending on the breaking 
pattern. We expect therefore to be the last chiral 
order parameter to vanish, after or together with all 
the other parameters (e.g., E), and thus 1^(3) < 1. The 
analysis of some properties of the small Dirac eigenval- 
ues suggests the same conclusion 
Even though this theoretical expectation has not been 
checked experimentally yet, we dismiss the case Y{3) > 
1 (unphysical in our opinion) in the remainder of this 
article. 

On the plots, very small values of X{3) do not neces- 
sarily correspond to almost vanishing pion masses. In- 
deed, we take into acount at least two different sources 
of chiral symmetry breaking : quark condensation IJ{3) 
and vacuum fluctuations L4, Lg. For X(3) = 0, chiral 
symmetry breaking may still be triggered by the fluc- 
tuations encoded in Lg, and the simulated pion mass 
may remain rather massive for m ^ . 
The pseudoscalar masses Mp are not the best cairdi- 
dates to observe a curvature due to chiral logarithms. 
The m-dependence is hardly different from a polyno- 
mial one. In addition, X(3) must be much smaller than 
1 to observe a variation in the curvature. 



Mi 



Ml . ...p 



—y- 

2L3 ^ Ljp ' 



where: 



L 



P ■ 



We obtain for the decay constants: 



- + 2]Ll + 8Ll 



(34) 



(35) 



(36) 



1 



4F2L3 



[4ct^ + 2dK] 



P2 



(37) 



A/2 



1 



4F2L3 



3. 3. 

-CT^ + 6(Jk + 2'''») 



where: 



CTp 



E 



t xiP + Ml 



In a similar way, we obtain for the masses: 



FX 



= gr|X(3) 



(38) 



(39) 



16 ( - + 2 ) 



16L^ 



r(3) 



AF^L^ 

FlMl ^ qr fl 
F?M3 



M"^ 

-^cir[Y{i)Y 



3a-7r + 2aK + -f?,, 




X{3) 



16\- + 2]Ll 



(40) 




Fig. 1. F^, F^AI^ and (respectively upper, middle and lower plots) normalized to their physical values as functions 
of g = rh/ms. The first (second) column corresponds to Z{3) = 0.8 (0.4). On each plot, full, long-dashed and dashed curves 
correspond respectively to X(3) — 0.8, 0.4, 0.2. Thick (thin) lines are drawn for r = 30 (20). All NNLO remainders are neglected. 



y(3) 

4F2L3 



5, 



2.2 Bare expansions at finite volume 

The above NLO chiral expansions are not bare expansions 
yet. We must reexpand the physical masses in powers of 
quark masses wherever the dependence on the latter is 
polynomial. Therefore, we have to identify the nonana- 
lytic pieces in the tadpole term a, which are a logarithm 
due to ultraviolet divergences and a pole due to infrared 
divergences. 



Since the (finite- volume) tadpole sum and the (infinite- 
volume) integral share the same ultraviolet divergences, it 
is convenient to introduce their difference |2U| : 



-y 



1 



(£2 + A'py 



(41) 



1/2) 



^ can be evaluated as an integral of known mathematical 
functions |19| : 



1 



(47r)3/2r(£ 



(42) 




Fig. 2. F^, F](Mx and (respectively upper, middle and lower plots) normalised to their physical values as functions 
of g = fh/ms. The first (second) column corresponds to Z{?,) = 0.8 (0.4). On each plot, full, long-dashed and dashed curves 
correspond respectively to X(3) = 0.8, 0.4, 0.2. Thick (thin) lines are drawn for r = 30 (20). All NNLO remainders are neglected. 



In order to turn the previous NLO chiral expansions 
into bare expansions, we can isolate the nonanalytic pieces 
in dp and expand the rest in powers of quark masses: 



dp _ Ml Ml 

L0(M2) £p 



+ ^1/2 (i,Mp) 



1 



LO(M|,) 
87r2 



log 



L^Mp 
LO(Mj,) 
Ml 



Cu2iL,L0{Ml)) 



l^Jlo(mI' 



where LO[Af|>] denotes the leading-order contribution to 
the meson mass: 





L0[M2] 


= qrY{3)M^ , 


(45) 


(43) 


LO[M|.] 




(46) 


(44) 


L0[M2] 




(47) 



In Adp. JXI we check that all the nonanalytic depen- 
dence on Ml in dp comes from the first two terms in 
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eq. H44|l . whereas the last (bracketed) is analytic in LO(Afp). 
We could take eq. H44|l as the bare expansion of dp. How- 
ever, a further simplification can be performed : the (last) 
logarithmic term in Eq. H44|l is small whatever the size of 
the fluctuations and contributes to ctp only at NLO. We 
choose to absorb these terms in the NNLO remainders of 
the masses and decay constants and to define the bare 
expansion of irp through: 



crp 
L3 



lo(m; 



87r2 



+ 



log 



1 



(48) 



ei/2(L,LO(M2)) 



L^JLOiAf 



where the ellipsis denotes 0(p*) terms that are absorbed 
in the NNLO remainders. 

First, let us notice that we recover the results of the 
previous section in the large- volume limit: only the first 
term in eq. 148|l survives, which is exactly the bare expan- 
sion of the infinite-volume tadpoles discussed in Sec. 11.11 
In addition, this choice for the bare expansion settles re- 
lated issues concerning the convergence of the expansion 
at finite volume. xPT in a finite box can be consistently 
formulated in three different regimes, called p-, e- and S- 
expansions, depending on the relative sizes of the box sides 
and the pion mass 17 . In particular, the p-expansion used 
here holds if: 



tions to the infinite-volume estimates: 



A{X) 



X{L) - X(oo) 

X{oo) 



(50) 



We obtain the following results for the decay constants: 



1 



1 



2S 



K 



(51) 



4F|l 



oo) 



3^ ■ 
2"" 



where: 



ei/2(i,L0[Af|.]) 



L3 



LO{Ml] 



L^Mp 



In a similar way, we obtain for the masses: 



A{FlMl) 



2S 



4F2(oo)M2(oo) 

X 3^, 

(g + l)ry(3) 

8F|.(oo) Mi(oo) 



K 



Ml 



1 ^ 



(52) 



(53) 



(54) 



3^ 
2 



-I- i^K + -z^n 



i^^L > 1 , M^L > 1 . (49) 

If the second condition is violated, the pion is too large to 
be contained in the box : the breakdown of the p-expansion 
is flagged by infrared divergences. In the case of small fluc- 
tuations, quark condensation is responsible for Nf = 3 chi- 
ral symmetry breaking. The second condition in eq. H49|l 
is translated into 2mBoi^ = LO(M^) x > 1. If cither 
m or _Bo are too small compared to the size of the box, 
the condition is violated and infrared divergences occur in 
dp. 

In the case of large fluctuations, the second condition 
in eq. 149|l has a different interpretation. We do not as- 
sume any specific mechanism for Nf — 3 chiral symme- 
try breaking. For non vanishing quark masses, Bq may be 
small without leading to a small pion mass : for instance, 
in the limit case i?o ^ 0, a significant pion mass can be 
generated by the vacuum fluctuations encoded in Lg and 
L4. Therefore, a small pion mass is obtained only if the 
quarks are light enough, and the p-expansion must break 
down for m too small but not for Bq too small. One can 
check that the proposed bare expansion of dp exhibits 
an infrared divergence when m ^ 0, but not i?o — > 0. In 
other words, even a very small three-flavour quark conden- 
sate does not lead to the breakdown of the p-expansion if 
vacuum fluctuations are significant enough. 

We could now take the NLO finite- volume expressions 
for the masses and decay constants, perform the replace- 
ment in eq. and study the resulting expressions. How- 
ever, we are mainly interested in the finite-volume correc- 



2.3 Size of the finite-volume corrections 

We begin with a sample plot of the finite-volume correc- 
tions. In fig. 01 A{Fl) and A{FlM'l) are plotted as func- 
tions of q for the specific choice of parameters r = 25 and 
Z(3) — 0.8. The left (right) column corresponds to L = 1.5 
fm (2.5 fm), corresponding to M^L = 1.1 {M^L = 1.8). 
The thick full, full, long-dashed and dashed curves corre- 
spond to X{3) = 0, 0.2, 0.4, 0.8 respectively. For X(3) = 0, 
the finite- volume correction to F^M^ vanishes at this or- 
der, see eq. 

We see that the corrections are very significant for 
L — 1.5 fm, but much smaller for L = 2.5 fm. In ad- 
dition the corrections are smaller for F^M^ than for the 
decay constant. Indeed, the finite- volume chiral expansion 
oiF^M^ contain an additional factor LO(M^) multiplying 
the tadpole terms dp. This tames the infrared divergences 
occurring when the pion mass is too small compared to 
size of the box. This damping factor is absent in the ex- 
pansion of the decay constants. 

We want the finite- volume corrections to remain small 
over a reasonable range of variation for q. Noticing that 
the size of the corrections is weakly dependent on Z{3), 
we introduce the quantity: 

DiX) ^ Max \AiX)\ | ^^^3 |(^)-<\ , (55) 

whose size we investigate as a function of r and 1^(3) for 
L = 2 and 2.5 fm, corresponding respectively to M^^L = 
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Fig. 3. The relative finite- volume correction to (upper plots) and F^M^ (lower plots). The left (right) column corresponds 
to L = 1.5 fm (2.5 fm). We set r — 25, Z{3) — 0.8. The thick full, full, long-dashed and dashed curves correspond to X(3) = 
0, 0.2, 0.4, 0.8 respectively. 



1.4 and Mt^L — 1.8. In figs. 01 and black regions corre- 
spond to D smaller than 5%, increasingly lighter regions 
to D smaller than 10,20 and 40% respectively. 

We observe the expected decrease of the finite- volume 
effects when the size of the box increases. Once again, we 
notice that the corrections to FpMp are much smaller 
than those to the decay constants, due to a better be- 
haviour at the approach of the infrared region. FpMp 
[P — TT, K) in large volumes {L = 2.5 fm) is a quantity 
for which we manage a good control of finite-volume ef- 
fects. 



we estimate the size of NNLO remainders in the physical 
case by attributing a 30% (10%) effect to an rus (m) fac- 
tor, which leads to 0{m'^) ^ 10%, 0{mms) ^ 3% and 
0(m2) ~ 1%. 

For lattice simulations where m varies, we take the 
following estimate: 



CK = 0{rhms) - O.lOg^/^ 
= O(m^) - 0.10, 



(56) 
(57) 



3 Constraining three-flavour order parameters 
on the lattice 

For definiteness, we take a volume of size L — 2.5 fm 
and set r — 25. However, we must keep in mind that 
the latter is a parameter which may vary between 15 and 
approximately 35. The qualitative conclusions that we are 
about to draw do not depend on the exact values of r, but 
one observes small shifts in the following plots when r is 
varied in its range. 

We take into account the remainders associated with 
the simulated masses and decay constants ep and dp (but 
not the indirect remainders d, d', e, e'). As a rule of thumb TS' 



in order to recover the physical case, i.e. 0{mms) ~ 0.10 
for m = TOs {q = 1) and O^rhiris) ~ 0.03 for rh = m 
(q^l/r). 

Figs. and [7| show the pseudoscalar decay constants 
and masses as functions of the quark mass ratio q — 
rh/nis- The left column corresponds to Z{3) = 0.8, the 
right one to Z{3) = 0.4. On each plot, the various bands 
show the impact of NNLO uncertainties for different val- 
ues of X{3) [fuU: 0.8, long-dashed: 0.4, dashed: 0.2]. We 
plot the results only when the finite-volume effects are 
smaller than 10%, but we do not include these corrections 
in the quantities plotted. 

At small q, the pion mass becomes too light, the finite- 
volume effects become larger than 10%, and we cannot 
trust the corresponding results because of potentially large 
^igher-order corrections (in such a case, we set the result 
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Fig. 4. The maximal finite-volume corrections D{F^) (upper plots) and D{F^M^) (lower plots) as functions of r and Y{3). 
The left (right) column corresponds to L = 2 fm (2.5 fm). Black domains correspond to D smaller than 5%, increasingly lighter 
domains to D smaller than 10,20 and 40% respectively. 



to on the plots). In particular the upper left plot in 
fig- El [F^/F^ for Z{3) = 0.8], where q must be rather 
large to tame finite- volume effects, confirms that the chiral 
expansions of the decay constants may suffer from sizeable 
uncertainties at small simulated masses because of large 
finite- volume corrections. 

From the previous analysis, the sensitivity to the three- 
flavour order parameters through the curvature seems more 



important in the case of the masses (which have smaller 
finite- volume corrections as well). We can isolate this ef- 
fect by considering the dimensionless ratios: 



Rtt — 



1 F^Mi 



Rk 



^'•^ (58) 



Figs. IHl and 1^1 indicate that these two ratios can provide 
a way of constraining the size of vacuum fluctuations by 




Fig. 5. The maximal finite-volume corrections D{F^) (upper plots) and D(F'^M]() (lower plots) as functions of r and Y{3). 
The loft (right) column corresponds to L — 2 fm (2.5 fm). Black domains correspond to D smaller than 5%, increasingly lighter 
domains to D smaller than 10,20 and 40% respectively. 



varying q from 1/r to 1. The larger r, the easier the dis- 
tinction between small and large fluctuations, even though 
the uncertainties due to NNLO remainders increase at the 
same time. 



4 Conclusion 

The presence of massive ss-pairs in QCD vacuum may in- 
duce significant differences in the pattern of chiral symme- 
try breaking between Nj = 2 and Nf = 3 chiral limits, i.e., 
when TOg remains at its physical mass and when m., is set 
to 0. This effect, related to the violation of the Zweig rule 
in the scalar sector, may destabilise three-flavour chiral 
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Fig. 6. and F^M^ (normalised to their physical value) as functions of q. The first (second) column corresponds to Z{3) = 0.8 
(0.4). On each plot, full, long-dashed and dashed bands correspond respectively to X{3) = 0.8,0.4,0.2. Cases where Y{3) > 1 
are not shown. 



expansions numerically, by damping leading-order (LO) 
terms proportional to 1^(3) and i^^(3), and by enhancing 
next-to-leading-order (NNLO) terms containing the O(p^) 
Zweig-rule violating low-energy constants L4 and Lg. In 
such case, a more careful treatement of chiral expansions 
is required to avoid uncontrolled corrections from higher 
chiral orders. 

In a previous work [ISj , we proposed a consistent frame- 
work to take into account the possibility of large vacuum 
fluctuations. Indirect hints from dispersive estimates 
I11II12| suggest that this effect might be significant, but an 
experimental determination has not been achieved yet. In 
this paper, we proposed to probe the size of ss-pairs fluc- 
tuations through lattice simulations with three dynamical 
flavours, with a strange quark at its physical mass nis and 
two lighter flavours of mass m. We focused on the masses 
and decay constants of the pions and kaons, and worked 
out chiral expansions which should exhibit small NNLO 
remainders even when LO and NLO terms compete nu- 
merically. The dependence of these observables on m can 
provide useful constraints on the structure of the Nf = 3 
chiral vacuum, and in particular on the size of vacuum 
fluctuations. Conversely, this dependence on vacuum fluc- 



tuations should stand as a warning about three-flavour 
chiral extrapolations on the lattice, which might prove 
more delicate to handle than usually assumed if vacuum 
fluctuations of ss pairs are signiflcant. 

We have also estimated the corrections due to the fi- 
nite spatial dimensions used in the lattice simulation. As 
expected, large volumes (L 2.5 fm) are required to pre- 
vent these corrections from spoiling the predictivity of the 
chiral expansions in the p-regime. These corrections do not 
seem enhanced in the case of large vacuum fluctuations for 
the observables considered here. 

Finally, we isolated two dimensionless ratios based on 
F^M^ and F^M]^ showing interesting features. They are 
not affected very strongly by finite-volume corrections, 
and their dependence on m could provide interesting in- 
sights at the size of vacuum fiuctuations. Therefore, it 
would be rather interesting to perform a lattice study of 
these ratios with three dynamical flavours, a reasonably 
large spatial box, and an action with good chiral prop- 
erties. The results might shed some light on the pattern 
of Nf = 3 chiral symmetry breaking and the low-energy 
dynamics of QCD. 
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Fig. 7. and F^Mj^ (normalised to their physical value) as functions of q. The first (second) column corresponds to Z{3) = 0.8 
(0.4). On each plot, full, long-dashed and dashed bands correspond respectively to X{3) = 0.8,0.4,0.2. Cases where Y{S) > 1 
are not shown. 
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A Bare expansion of finite-volume tadpoles 

The finite- volume tadpole ap is related to its infinite- 
volume counterpart through ^1/2 which can be reexpressed 
as: 

Ci/2(i, Af2) = 2 hm g['^ = ^gfA > (59) 

where gi"^^ has been introduced in App. A of ref. 20 , and 
the zero-temperature limit is given in eq. (B.l) of the same 
reference ^. 



^ We take this opportunity to correct a typo in this equation, 
which should read: 



(d-i) 

,. (d) 9,-1/2 

hm g;. ' = — 



(60) 



From App. B of ref. [201 j it is straightforward to de- 



(3) 

termine the expansion of g\ L in powers of M: 



yi/2 



7n 
71 



47ri^ 



1 

47r 

1 

47r 



27r Af2L2 



2t2\ 



ML 47r 



loglM-'i 



(61) 



n=0 

an+1/2 + 



(2n+ l)(?i- 1) 
5 



^3/2 - r'(l) ~ log(47r) - 



(63) 



with q;„_|_i/2 are known numerical coefficients (the first few 
are displayed in table 1 of ref. (201 )■ 

This leads to the following expansion of ^1/2 (i, M): 



6/2(i,M2) 



M2 



ML3 



57r 

00 



log(M-'L^) 



(64) 



27r ^ n! 




Fig. 8. The ratio R-n as a function of q. The first (second) column corresponds to r = 20 (30). The first (second) row deals 
with Z{3) = 0.8 (0.4). On each plot, full, long-dashed and dashed bands correspond respectively to X{3) — 0.8,0.4,0.2. Cases 
where Y{3) > 1 are not shown. 



Thus, the nonanalytic dependence of (jp on LO (Af|>) conies 
only from the pole and the logarithm singled out in eq. (|44(l . 
According to our prescription, the bare expansion is ob- 
tained once the polynomial terms are reexpanded in pow- 
ers of quark masses. At our order of accuracy, this amounts 
to replacing Mj, by LO(Af|>) in the polynomial pieces, i.e.: 



LO(M|.; 
87r2 



log[LO(Af|,)L2 



Ci/2(2.,LO(A?2))_ 



L0(AfpL3) 



X3 



Ml Ml , 
^log^+Ci/2(i,M^) 

L0(Af2) , Ml 
log 



87r2 



[ MpL3 



LO(A^|,) 



87r2 

hO{Mj 



\og{MlL^) 



Stt' 



log^ 



M| 
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MpL^ 



LOjM^) 



log(Af|,L2 



LO(A? 



(65) 
(66) 



Mj 
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LO(Afp) LO{Ml) 



87r2 



■log- 



A//2 



ei/2(i,L0(Af|,)) 



(67) 
(68) 



(69) 



L0(Afp)L3, 



The bracketed term has only a polynomial dependence 
on L0{MD. This expression yields eq. (|48|l and the bare 
expansion of crp. 




Fig. 9. The ratio Rk as a function of q. The first (second) column corresponds to r = 20 (30). The first (second) row deals 
with Z{3) — 0.8 (0.4). On each plot, full, long-dashed and dashed bands correspond respectively to X{3) = 0.8,0.4,0.2. Cases 
where Y{3) > 1 are not shown. 
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